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Abstract 



The relation between the pole quark mass and the MS'-renormalized mass is governed 
by an infrared renormalon singularity, which leads to an ambiguity of order Aqc d in the 
definition of the pole mass. We use the renormalization group and heavy quark effective 
theory to determine the exact nature of this singularity up to an overall normalization. 
In the abelian gauge theory, the normalization is computed partially to next-to-leading 
order in the flavour expansion. 
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1. It is well known that perturbative series in QCD diverge. Recently there has been 
a resurge of interest in learning about nonperturbative effects in a variety of situations 
through this divergence. Among the quantities that have become an object of scrutiny is 
the pole mass of a quark, defined perturbatively as the location of the singularity of the 
renormalized quark propagator. The pole mass can be related to the MS renormalized 
mass - which in principle can be measured to any accuracy at a very high energy scale 
- by the seriesQ 



oo 

™pole = TBj^g + r n& n+1 (1) 
n=0 

where the first two coefficients r$ and r\ are known Jl|. Being infrared (IR) finite and 
scheme- independent [[|, the pole mass is the natural mass parameter in heavy quark 
systems, whose scale is governed by the mass of the heavy quark. Moreover, the pole 
mass is implicit in the construction of heavy quark effective theory (HQET), which has 
become the major theoretical framework for the description of heavy hadron decays ||. 

Recently, investigation of the set of diagrams symbolized by Fig. la revealed that 
in large orders, the coefficients r n are strongly contributed by IR momenta in the 
integration over the gluon line and diverge as j|, [5| 

rn ^^ e 5/a M( _ 2A) n n , (2) 

where (3q is the first coefficient of the /3-function [We use (3(a) = ji 2 {da) / \dji 2 ) = 
p a 2 (l + En=i c n a n ). C F = 4/3 for SU(3) and C F = 1 for U(l).]. If one attempts to 
sum the divergent series with the help of the Borel transform, the Borel sum 



"ipole = ^MS + / dt B [ m po\e](t) B [mpole] (*) = Y r n — , (3) 



is ambiguous of order Aqcd from a simple IR renormalon pole at t — — 1/ (2/3 ) - In the 
following we adopt the convention that all quantities that have divergent expansions are 
understood as Borel sums as in Eq. with the contour deformed above the singularities 
on the positive axis. Then the ambiguity of order Aqcd is equivalently expressed as an 
imaginary part of m po i e of this order. 

The ambiguity in the definition of the pole mass can immediately be translated 
into an ambiguity of the parameter A# Q , which appears in HQET, if defined as usual 
through the expansion of the heavy hadron mass 

m HQ = m po i c + A Hq + 0(l/m po i e ). (4) 

Physically, this ambiguity of Ah q arises, because the contribution of the spectator to 
the mass of the hadron can not be separated from self-energy corrections to the heavy 
quark, once these become sensitive to the long range part of the Coulomb potential. 



1 A normalization point [i is understood for m-rjrs and the MS coupling a 
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This fact is known in a different language to the lattice community, where the extraction 
of the constant term in the heavy quark potential has faced difficulties precisely because 
of the contamination from self-energy contributions. 

In this letter, we address a more formal issue and explore to what extent the large 
order behaviour of the coefficients r n can be determined without any approximation. 
In general, one expects the behaviour of Eq. (Q) to be modified to 

r n n =°° N v (-2/?o)" T(n + 1 + b) ^ + £ , (5) 

turning the simple pole of B[m po \ e ](t) at t = — l/(2/?o) into a branch point of strength 
1 + 6. We shall determine the exact nature of the singularity, i.e. b and the Sk, entirely 
in terms of the /3-function. The analoguous coefficient b is known to date only for the 
QCD corrections to e + e~ annihilation and the GLS sum rule || [7]], where, however, the 
corresponding singularity is not responsible for the dominant divergent behaviour of the 
series. In these cases one may exploit the short distance expansion of the appropriate 
current product and a relation to higher dimension operators. The short distance 
expansion is not available for the pole mass. Instead we will utilize the observation of 
Ref. f|], that the IR renormalon at t = — l/(2/3o) in the pole mass is closely related 
to an ultraviolet (UV) renormalon singularity at the same position in the self-energy 
of a static quark. The normalization iV remains elusive. We calculate it partially to 
next-to- leading order in the flavour expansion in the abelian theory. 

Though we doubt that our results can be used phenomenologically to improve the 
perturbative relation Eq. (P, we believe that the arguments that lead to the charac- 
terization of Eq. (^) are interesting by themselves and may shed more light on the 
simplicities and complexities of divergent series in renormalizable theories. 

2. In this section we consider the self-energy of a static quark, described by the 
leading order effective Lagrangian of HQET. Its Lorentz covariant form is 

£ c s = h v iv ■ Dh v + Aight, ( 6 ) 

where v is the quark velocity and £u g ht the QCD Lagrangian for the light degrees 
of freedom. In contrast to the self-energy of a quark of finite mass, the self-energy 
of a static quark is linearly UV divergent. Starting from a bare Lagrangian of the 
above form, renormalization will in general involve a counterterm h v h v , such that the 
renormalized Lagrangian contains a residual mass term 5m res oc [i for the static quark. 
We choose a minimal subtraction (MS) like subtraction scheme (i.e. MS or schemes that 
differ from MS by a change of scale /i). Then a residual mass does not arise to all orders 
in perturbation theory, but the subtractions are such that the series expansion of E e g 
has an UV renormalon divergence analogous to Eq. (|5|) H . We emphasize the choice of 
MS-like schemes. Large order coefficients are scheme-dependent and in general (mass- 
dependent) schemes, this particular UV renormalon divergence can be eliminated at 
the price of introducing a residual mass, which is adjusted precisely to this divergence. 
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At this point we return to the old observation by Parisi M that the imaginary parts 
of the Borel-summed Green functions due to an UV renormalon at position t = n/ (3 in 
the Borel plane are proportional to a sum over zero-momentum insertions of all possible 
local operators of dimension In + 4, which can be considered as a simple version of a 
factorization theorem. Generically 

ImG(p fc ) = ^ Yl EdiGoAPk)- (7) 

d AimOi=d 

Although we are not aware of a general proof of this statement, we note that it has 
been confirmed in a nontrivial case of the photon vacuum polarization || [K| as far 



as dimension six operators are concerned. In fact, Eq. (|?p is almost intuitive, if one 
understands the emergence of UV renormalons in the context of renormalization in 
general. After the theory is regulated with a dimensionful cutoff, renormalization can be 
expressed as hiding all divergent terms in a large-cutoff expansion into a redefinition of 
the low-energy parameters. The remaining cutoff-dependence determines the remaining 
sensitivity to small distances, which is suppressed by a power of the cutoff. It is these 
"left-overs", which give rise to divergent series, when the series is expressed in terms 
of the renormalized low-energy coupling. However, in principle, order by order in the 
inverse cutoff the remaining sensitivity to small distances could be removed by adding 
local higher dimension operators to the Lagrangian, a procedure that is indeed widely 
used in the construction of improved lattice actions. By the same token, the imaginary 
parts created by UV renormalons can be described by Eq. (|7|). 

In the particular case of the self-energy, the situation is in fact much simpler. The 
UV renormalon at t = — l/(2/?o) is related to the dominant (linear) cutoff-dependent 
term and one does not have to deal with a small left-over from renormalization. As 
indicated by the position of the pole and the discussion above, the relevant operators 
have dimension three. There exists only one dimension three operator, h v h v , and the 
general expression Eq. (0) reduces to 



Im See (vk, fx, a) = E(n,a) S^ ^^ ' a ) ' ( 8 ) 



as far as the renormalon at t = — 1/(2/3 ) alone is concerned. Here ^ e s,h v h v denotes the 
self-energy with one zero-momentum insertion of h v h v . Eq. (§) is illustrated to leading 
order in Fig. 2. Comparing the renormalization group equations for £ e ff (and therefore 
its imaginary part) and ^ e Sh v hvi one finds that the coefficient satisfies 



2 d . „, s d 



The anomalous dimension 7^,^(0:) of the operator h v h v vanishes to all orders in per- 
turbation theory, because the operator is related to the conserved current h^^hy by 
spin symmetry. A different way to see this is to observe that mqh v h v is renormalization 
group invariant, but the HQET expansion parameter tuq does not run. By dimensional 
arguments, the explicit /i-dependence of E must be an overall factor /1 and the solution 
of Eq. (SI) is 
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E (p, a) = const x fi exp I y da' ^g^rj 1 = const x Aqcd- (10) 

Together with Eq. (§) this determines completely the UV renormalon singularity of the 
self-energy up to an overall constant in terms of the /^-function and the perturbative 
expansion of S eff Note that it depends onvk/p only through 1 /n-corrections to the 
leading asymptotic behaviour in accordance with the expectation that the divergence 
of the derivative of the self-energy is suppressed, because it is linearly divergent only 
through linearly divergent subdiagrams. For later use, we note the equality 



which can easily be shown diagrammatically by routing the external momentum k 
through the heavy quark line. The unity on the r.h.s. appears, because the tree dia- 
gram on the r.h.s. of Fig. 2 is included in the definition of S eff ^^. 

3. We now relate the UV renormalon in the self-energy of a static quark discussed 
above to the leading IR renormalon in the pole mass, generalizing the leading order 
observation of Ref. M. To this end we sandwich the full inverse propagator S* -1 of a 
quark in QCD between two projectors P v = (l+^)/2 and define 

P v S' 1 (p,m^g,fJ,;a) P v = P v Sp l (vk,m Q , (j,;a) (12) 



with p = rriqv + k. As long as p is real and away from the zero of S" _1 , (the Borel sum 
of) S~ l does not have an imaginary part of order Aqcd, because its Borel transform 
has no singularity at t — — 1 / (2/5 ) • Odd powers of Aqcd can arise only from on-shell 
Feynman integrals or expansion around mass-shell. For the purpose of our argument, 
we shall assume the hierarchy mq 3> vk 3> Aqcd- Then the largest imaginary part of 
^p 1 comes from a singularity of its Borel transform at t — —1//3q, which is or order 
AQ CD /vk, parametrically smaller than Aqcd- 

We wish to expand Sp -1 in vk/rriQ about vk = 0, where p, mq and vk should be 
real. Anticipating an imaginary part of the Borel sum of m po i e of order Aqcd, we define 

rriQ = m po i e — 8m 5m = iSm = iImm po i e . (13) 



To all orders in perturbation theory, Sp 1 can be matched onto the self-energy of a static 
quark computed from HQET without a residual mass term, Eq. @, and expansion 
parameter rriQ = m po \ c . Thus defining k! = k — 5mv and using 5m <C vk, one obtains 



Sp 1 (vk,m Q ,fi) 



C 
C 



K-Eeff K,/x))+o 

(?) 



~'|2 



^pole 



vk — E e g (vk, fx) — 5m I 1 



as 



off 



dvk 



(vk/n) 



(14) 
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+ *™C'(m) [vk-Y> cS (vk^)} + 

n \ n J 

with a matching coefficient C that contains all mg-dependence. While E e g- contains an 
imaginary part of order Aq C d, the l.h.s. of Eq. QI3j ) does not. Investigating the scaling 
of all potential imaginary parts in quantities on the r.h.s. of Eq. ([TJ]) with vk, // and 
Aqcd, one finds that the condition that imaginary parts of order Aqcd cancel, can be 
expressed as 



Im 



vk - S eff (vk, fi) -Sm(l- ^^-(vk/fi) 



0. (15) 
Combining this condition with Eqs. @ and (|TTD , we obtain the final result 

Imm po ie = 5m = —E((j,,a) = const x Aqcd (16) 

2 ~ c 
2/?o 



const x ^e 1 /(2/3o«) a d/(2A,) ( 1 + c J^A a + . . 



where E is defined by Eq. (^) and given by Eq. (PD[). The second line of Eq. ( |To| ) is 
easily translated into the form of the singularity of B[m po ie](t) at t = — 1/(2/3q), and 
therefore into the large order behaviour of the coefficients r n , Eq. (f|), which relates the 
pole mass to the MS mass^]. From Eq. (pTOj ) or ( Jl6| ) one deduces that except for the 
overall normalization iV all quantities that characterize Eq. @ are determined by the 
/5-function. For instance, 

6 = -2A S1 = -ST' (17> 

and, in general, l/ra fc -corrections to the leading asymptotic behaviour involve the f3- 
function coefficients up to Ck+i- We stress that these conclusions follow with almost no 
dynamical input. We used essentially only the linear UV divergence of the self-energy 
of the static quark. The only "nonperturbative" parameter that remains and is not 



fixed by very general arguments, is the overall normalization [|Tl], |12 |. This is because 
a nonperturbative approach comparable to the semiclassical expansion for instanton- 
induced (combinatorically) divergent series is lacking for renormalons. 

For the ambiguity present from the perturbative definition of the pole mass, Eq. ( Jl6| ) 
shows that it is proportional to Aqcd with no multiplying logarithm and series in 
ln(/i / Aqcd) > which is in fact the only possible form compatible with scheme-invariance of 
the pole mass. In our formal derivation, the absence of logarithms is due to the vanishing 
anomalous dimension of h v h v , which indeed is directly related to the renormalization 
group invariance of the pole mass as noted above. The bonus of the formal derivation 



2 For this argument, we assume that no "unknown" singularities are present. Then the IR renor- 
malon at t = — l/(2/3o) lies closest to the origin of the Borel plane and determines the asymptotic 
behaviour of the perturbative coefficients. 
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is Eq. (0), which allows to calculate (approximately) the overall normalization from the 
self-energy in HQET. 

As mentioned before, the renormalon ambiguity in m po i e translates into an ambigu- 
ity in the A-parameter of HQET, see Eq. (|j). It is tempting to estimate this ambiguity 
from a leading order calculation of the normalization. The arbitrariness of this proce- 
dure is reflected in the different values, spanning <5m po i e = SA Hq = (50 — 300) MeV, 
quoted in Refs. [|, |5| and the present considerations can not improve the expectation 
that this ambiguity is of order Aq CD . 

4. In a nonabelian theory, a systematic approach to calculate the normalization of 
renormalons is not known. Explicit calculations can be performed in an abelian toy 
model (QED with Nf massless fermion for Aight) i n an expansion in the number of 
flavours^. To organize this expansion, one defines a = Nfat and expands in a and 1/Nf. 
The /3-function for the rescaled coupling is f3(a) = 6o fl2 (l + c±a + ...), 6 = 1/(3tt) and 
Ci = 3/(47r). We also expand the normalization 

where Nq is minus the asymptotic behaviour of coefficients obtained from the set of 
self-energy diagrams in Fig. la. To next-to- leading order in 1/Nf, one has to extract 
the asymptotic behaviour from the diagrams symbolized by Fig. lb-d, which according 
to Eqs. @ and (K7h is of the form 



1 + — { — — In n + Ni 

Nf\ 2(3 




x(l +0 (!)). (i c„ 



A simplification arises, because diagrams lc and Id do not contribute to the overall 
normalization, when only one of the two gauge boson momenta is large. For example, 
diagram lc amounts to a radiative correction of the r.h.s. of Fig. 2 and therefore affects 
the asymptotic behaviour of £ e // only at the level of 1/n-corrections, whereas diagram 
Id does not contribute to the UV renormalon at t — — 1/(2/3q) at all, because it has 
no linearly divergent subdiagrams. The present case is markedly distinct from - and 
simpler than - the vacuum polarization in QED and its first UV renormalon at t = 1/(3q, 
where diagrams with two chains of fermion bubbles and one large internal momentum 
are not suppressed ||. The reason is that for the self-energy, we are interested in the 
dominant large-momentum contribution to the Feynman integrals. A contribution to 
N± from diagrams lc and d can only come from the region, where both gauge boson 
momenta are large and of the same order. 

Diagrams la and b are given by the first two terms in the 1 /Ay-expansion of diagram 
la with integration over the full gauge boson propagator. The Borel transform (defined 
analogous to Eq. (Q)) of the class of diagrams obtained by insertion of the full propagator 
is 



3 The same expansion can also be applied to the nonabelian theory, but can not be taken seriously 
numerically, since it is singular at Nf — 33/2, when (3q is zero. 
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AttCf f d 4 g (vk) 2 -q 2 



iN f J (2vr) 4 v(k + q) q' 



B 



a 



l + U{-q 2 /fi 2 ,a) 



(t), 



(20) 



as discussed in a general context in Ref . JTT| . Expanding the gauge vacuum polarization 



n. 



U(q 2 )=U (q 2 ) + l-U 1 ( q 2 ) + 



n (g 2 



-b$a\n.X X 



9 ' 



(21) 
(22) 



one finds 



B 



a 



l + Ii(q 2 ) 



(t) = X bot / ds sX boS B 

N fi 



ni 



O: 



(t~s) + ^ 



(23) 



where the first term corresponds to diagram la and the second to lb. The Borel 
transform of Hi/ a is known [13|] and can be written as 



B 



Hi 



a 



(t) = X bot F(t) - G(t), 



(24) 



where F(t) is a scheme-independent function given by [|13|, |14 

1 g (-i)k k 

k=2 



F(t) 



tfl-(l + t {) 2 £(l- 



i + M) 2 ) 2 " t' 



(25) 



The divergence of F(t) at t = is regulated by G(t), which is scheme-dependent 
and without poles in t. But, since a/(l + IT(g 2 )) is scheme-independent in QED, the 
scheme-dependence of its Borel transform is compensated by the scheme-dependence 
of the coupling in the exponent of the inverse Borel transform. Consequently, the 
subtraction dependence of G(t) can be expressed in terms of the /3-function. Define 
F reg (t) = F(t) + &i/t, G reg (t) = G(t) +ci/t, such that both F Teg (t) and G ieg (t) are finite 
at t = and 



Z(u) 



4ttC f f d 4 q (vk) 2 -q 2 



AttC f 



N 



vk 



f 



(2vr) 4 v(k + q) 
2vk^ ~ 2u 
I' 



x- 



(26) 



— ) e 5u/3 (-6) 



r(-i + 2tt)r(i 



u 



r(2 + u) 



Then 
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5[£ e a+b ](i) = /(-M) 



t t 
I - — J ds sF reg (t - s) + — J ds sI(-b s)G reg (t - s) 



t 



+ WJ ds^(A-M)-/(-M)- 

* o 



(27) 



It is simple to extract the leading singular behaviour as t — > — 1/ (26 ) from this equation. 
The term involving G reg produces a In (1 + 26q£) singularity and contributes only to 1/n- 
corrections in Eq. ([19|), consistent with the scheme-independence of N±. The result is 



t-+-l/(26o) C F 



vriV 



3 5/6 M. 



/ 



1 + 26 i 



1 + 



2b 



ln(l + 26 t) + ^| + 



(28) 



with 



-l/(26 ) 



ri /" 27 4 

A ^ 2T ~ J dssi? re g (-l/(26o)-s) = yln- 



(29) 



+12 E 

fc=2 



(-l) fc fc 
(fc 2 - l) 2 



1 

21 2 



In — — H ^5 — ^artanh A; 1 — artanh (2/c) 1 



fc 2 - 1/4 



2fc 3 



The last step consists in converting the singularity of the Borel transform in Eq. ( p8|) into 
the asymptotic behaviour of the coefficients r n . We reproduce correctly the coefficient 
of the Inn-term in Eq. (|T9D, which can already be deduced from Eq. (Pj) without any 
calculation, and find (7^ is Euler-Mascheroni's constant) 



26 



= 0.6831... 



(30) 



This concludes the partial evaluation of the normalization in next-to-leading order in 
the abelian theory. Note that to obtain N\ one needs the exact vacuum polarization III, 
in accordance with general expectations [11]. We repeat that a potential contribution 
from diagram c and d from the region where both internal momenta are large is still 
missing. The presented calculation of N\ should be considered as an illustration how 
the most complicated diagram b can be easily evaluated, using known results on the 
vacuum polarization. The Borel transform of diagrams c and d is given by (u = —bot, 
v = -60s) 




r(i + 2v) 



T(2v) 



r(-i + 2«)r(i - u)r(i -u + v) 

T(2 + v)T(2 + u-v) 
T(-l + 2u-2v) . (31) 
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However, before one can use this expression, one has to extend the renormalization of 
the Borel transform [|J to diagrams with two chains. 

5. The pole mass or, respectively, the self-energy of a static quark is probably the 
simplest quantity, for which the leading renormalon can be determined up to a global 
normalization. Its simplicity might also open a path to tackle the outstanding problem 
to identify which diagrams in addition to fermion bubbles are responsible for converting 
the abelian fi to its nonabelian value, which is expected on general grounds in Eq. (|5[). 
To some extent, this issue can be addressed at order 1/Nj in the flavour expansion of 
the nonabelian theory. The obvious replacement of one fermion bubble in the diagram 
la by a gluon or ghost bubble is not sufficient to produce the l/iV|-contribution from 
expanding (3q . It would be interesting to check that the two- loop diagram involving the 
three-gluon vertex dressed by fermion bubbles can make up for the missing piece. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408380vl 



(a) (b) 




Fig.l (a) Leading order diagram, (b) - (d) Diagrams to next-to-leading 
order in the flavour expansion in the abelian theory. Each diagram (a) - (d) 
represents the set of diagrams obtained by inserting any number of fermion 
loops in any of the photon (gluon) lines. 




Fig. 2 Eq. (8) to lowest order. The gluon line is dressed by a chain of bubbles. 




